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Abstract 

We describe how to compute the intersection of two Lucas sequences of the forms 
{U„(P, ±1)}^ =0 or {V„(P, ±1)}£L with P e Z that includes sequences of Fibonacci, Pell, 
Lucas, and Lucas-Pell numbers. We prove that such an intersection is finite except 
for the case U n (l, -1) and U n (3, 1) and the case of two V-sequences when the product 
of their discriminants is a perfect square. Moreover, the intersection in these cases 
also forms a Lucas sequence. Our approach relies on solving homogeneous quadratic 
Diophantine equations and Thue equations. In particular, we prove that 0, 1, 2, and 5 
are the only numbers that are both Fibonacci and Pell, and list similar results for many 
other pairs of Lucas sequences. We further extend our results to Lucas sequences with 
arbitrary initial terms. 

In 1962 Stein fTT| using elementary arguments proved that if two Fibonacci sequences 
(with different initial terms) share three terms then they are identical from some terms 
on. This result was later generalized by Revuz Q to arbitrary sequences satisfying the 
same second order linear recurrence. A similar problem is to determine the intersection 
of two sequences satisfying distinct linear recurrences. Mignotte, 1978 [6] proved that 
the intersection of two sequences is finite unless the roots of their characteristic polyno- 
mials are multiplicatively dependent. Matyas, 1981 [5j gave a criterion for determining 
whether two second order linear sequences have nonempty intersection. Posing the prob- 
lem as a system of Pellian equations, Pinch, 1988 |7j and Tzanakis, 2002 ITT511 proposed 
computational methods for solving such systems. 

In the current paper we focus on the computational matters of finding the intersec- 
tion and provide some explicit results for a class of Lucas sequences. While the general 
structure of such intersections is well known (6j [5]. [7|. [151, our approach (outlined in the 
proof of Theorem [6]) is different from previous studies (see 11151 for a review of various 
computational methods). In short, we reduce the problem to a finite number of Thue 
equations, each of which has a finite number of solutions. In practice the intersection of 



particular Lucas sequences can be computed with PARI/GP computer algebra system |13| 
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providing a functionality for solving Thue equations, based on Bilu and Hanrot's im- 
provement [lj of Tzanakis and de Weger's method ||T6| . We also characterize the cases 
when the intersection is infinite and show that it also forms a Lucas sequence. 

While the current paper was under review, our attention was drawn to the paper 
of Szalay, 2007 fl2|, which employs a similar approach for solving systems of Pellian 
equations. We therefore feel obliged to underline the differences between fl2| and our 
paper. In Theorem [6] we consider a general system of two homogeneous quadratic equa- 
tions in two indeterminants and formulate conditions under which such a system has a 
finite number of integer solutions. The systems of Pellian equations considered in |12| 
represent a particular case of the system from our Theorem [61 while the algorithm and 
the corresponding proof of solutions' finiteness in ||T2f are flawedQ In |[T2| it is used the 
result from [4] that the denominator in the general solution to a Diophantine equation 
Ax 2 + By 2 + Cz 2 = 0, derived from a particular solution (x , yo, z ) with z + 0, divides 
2A 2 BCZq. In Theorem [5] we prove a stronger result that this denominator actually divides 
2 lcm(A, B)CZq, leading to a smaller number of Thue equations to solve. This result also 
corrects an error in Corollary 6.3.8 of 0. We further focus on multiple issues specific to 
the intersections of Lucas sequences which are out of scope of |i~2f. Said that, we refer to 
fl2| for many examples of systems of Pellian equations and a review of their applications 
and related literature. 

The paper is organized as follows. In Section [T] we give definition and basic properties 
of Lucas sequences. We develop an algorithmic approach to solving systems of quadratic 
Diophantine equations in Section |2] and apply it to the problem of finding intersections 
of Lucas sequences in Section [3] Some numerical results are given and a number of open 
questions are posed in Section HI Finally, we discuss generalizations of our approach to 
Lucas sequences with arbitrary initial terms in Section [5] 



1 Basic properties of Lucas sequences 

The pair of Lucas sequences U„(P, Q) and V n (P, Q) are defined by the same linear recurrent 
relation with the coefficient P, Q e Z but different initial terms: 

U (P,Q) = 0, U 1 (P,Q) = 1, U n+1 (P,Q) = P-U n (P,Q)-Q-U n - 1 (P,Q),n>l; 
V Q {P,Q) = 2, V 1 (P,Q) = P, V n+l {P,Q) = P-V n {P,Q)-Q-V n . l {P,Q),n>l. 

Some Lucas sequences have their own names: 



Sequence 


Name 


Initial terms 


u»(i,-D 


Fibonacci numbers 


0,1,1,2,3,5,8,13,21,... 


Vn{l,~l) 


Lucas numbers 


2,1,3,4,7,11,18,29,... 


U n (2,-1) 


Pell numbers 


0,1,2,5,12,29,70,169,... 


V n (2,-l) 


Pell-Lucas numbers 


2,2,6, 14,34,82,198,... 



1 In [12] it is incorrectly claimed that the constructed quartic equation (2.5) represents a Thue equation. We 
found a counterexample with (a\, b\, C\) = (5, -1, -1) and (a 2 , b 2 , c 2 ) = (20, -1, 1) in (2.1) - (2.2) which yields 
(a,b,c) = (25,-1,-1) in (2.4). From a basic solution {Xq,Yq,Zq) = (1,0,5) it further leads to the reducible 
quartic polynomial T(s, r) = 3125s 4 - 2250sV + 5r* = 5(25s 2 - 20sr - r 2 )(25s 2 + 20sr - r 1 ) in the l.h.s. of (2.5). 



Other examples include Jacobsthal numbers U n (l, -2), Mersenne numbers U n (3, 2) etc. 
In the current paper we focus on the case ofQ = lorQ = -l. We also notice that 

U n (-P, Q) = (-l) n+l U n (P, Q) and V n (-P, Q) = (-l) n V„(P, Q) 

and restrict our attention to Lucas sequences with P > 0. Similarly, we exclude from 
consideration terms with negative indices, noticing that they may deviate only in signs: 

tL„(P, 1) = —U n (P, 1), V. n (P, 1) = V n (P, 1), 

tL„(P,-l) = {-lf +l U n {P,-l), V- n (P,-l) = {-l) n+l V n {P,l). 

From further consideration we also exclude the following degenerate cases: 



(P,Q) 


U n {P,Q) 


V M (P,Q) 


(0,1) 


u 2m = 


V 2m = 2 ■ (-1)'" 

Vjm+l = 


(o,-i) 


u 2m = o 

U2m+1 = (— l) m 


^ 2m = 2 

^m+l = 


(i,i) 


U 3 m = 
^3m+l = (-1)" 1 
L?3m+2 = ( _ 1)"' 


V 3m = 2 
V 3m+1 = (-If 

V 3m+2 = (-l) m+1 


(2,1) 


U w = m 


V m = 2 



It is easy to see that under the described restrictions both U n (P, Q) and V n (P, Q) are positive 
for n > 0. 

The characteristic polynomial of Lucas sequences {U n (P, Q)} and {V n (P, Q)} is A 2 -PA+Q 
with the discriminant D = P 2 - 4Q. For non-degenerate sequences, the discriminant D is 
a positive non-square integer. 

Let a = and jS = P ~ 2 V ^ be the roots of the characteristic polynomial, then the 

following explicit (Binet-type) formulas take place 

Q = «]6 

D = (a- ft) 2 

U n (P,Q) = ^ 

V n (P,Q) = a n +p n 

In particular, these formulas imply that 

V n (P,Q) 2 -D-lT„(P,Q) 2 = 4Q". (1) 

For |Q| = 1, it means that the pairs (V n (P, Q), U n (P, Q)) form solutions to the equation]^ 

x 2 - Dy 2 = ±4. (2) 

The converse statement can be used to prove that given positive integers belong to 
{V n (P, Q)} or {U n (P,Q)\ respectively (stated without a proof for Fibonacci numbers in 

my. 



2 Here and everywhere below ± in the r.h.s. of an equation means that we accept both signs as solutions. 



Theorem 1. Let P, Q be integers such that P > 0, |Q| = 1, (P, Q) gt (3, 1), and D = P 2 - 4Q > 0. 

If positive integers u and v are such that 

v 1 - Du 2 = ±4, 

then 

ue{U n (P,Q)} and ve{V n (P,Q)}. 

Proof. We notice that D is congruent to or 1 modulo 4, implying that the general solution 
to equation (0 has the form (x, y) = (Xk, yk), where 



x k + y k ^lD (xi + yi Vd 



and {x\, yi) is a solution with the smallest positive yi (see Proposition 6.3.16 in [2J). Since 
there is a solution (x, y) = (P, 1) with the smallest possible positive y = 1, we can take 

(xi,yi) = (P,l)B 

Therefore, for some «£Z 

u + u Vd /p+Vd\ 

— =± (— ) 

implying that 

W = ±LT n (P,Q) and v = ±V n (P,Q). 
Since u > and v > 0, we have 

M = U W (P,Q) and e = V|„|(P,Q). 

□ 



2 Quadratic Diophantine equations 

While theorems in this section can be viewed as consequences of the fundamental results 
of Thue |14| and Siegel IfTOl , we focus on computational aspects and give independent 
self-contained proofs of algorithmic nature. 

Theorem 2. A system of Diophantine equations 

fi{m,n) = k x 
f 2 (m, n) = k 2 

where f\(m,ri) and fi(rn,ri) are linearly independent homogeneous quadratic polynomials with 
integer coefficients and (ki,k 2 ) ^ (0,0) are integers, has a finite number of solutions. 

3 In the case (P,Q) = (3,1), there is also another solution (x, y) = (1,1) with y = 1, see Theorem [7] for 
further details. 



Proof. Without loss of generality assume that k\ + 0. 

If m orn is zero, then the other indeterminate may take only up to two values satisfying 
the quadratic equation fi(0, n) = ki or fi(m, 0) = k\. For the rest assume that both m and n 
are non-zero. 

The given system implies an equation g(m, n) = where g(m, n) = k 2 fi(m, n) - kif 2 (m, n) 
is a homogeneous quadratic polynomial. 

If g(m, n) = cmn for some coefficient c ± 0, then there are no non-zero solutions. 

If g(m, n) ^ cmn, without loss of generality we assume that the coefficient of m 2 in g(m, n) 
is non-zero, implying that g(m,n) = h(m/n)n 2 where h(x) is a quadratic polynomial. In 
this case, we have a quadratic equation h(m/n) = that determines up to two possible 
values of the ratio m/n. For a fixed value of m/n, the equation fi(m,n) = k\ has only a 
finite number of integer solutions. 

Therefore, the given system of Diophantine equations has only a finite number of 
solutions. □ 

Theorem 3. A Diophantine equation f(m, n) = k, where f(m, n) is a homogeneous polynomial of 
degree 4 with integer coefficients and k + Qis an integer, may have an infinite number of integer 
solutions only if f(m, n) = c • g(m, n) 2 , where c is an integer and g(m, n) is a polynomial of degree 
2 with integer coefficients. 

Proof. Since the polynomial f(m, n) is homogeneous, it can be represented as f(m, n) — 
n 4 g(m/n) where g(x) is a univariate polynomial with integer coefficients. Consider a fac- 
torization of g(x) over Q that implies a factorization f(m, n) into irreducible homogeneous 
polynomials over Z. 

If f(m, n) has two linear factors a^m + bin and a 2 m + b 2 n that are linearly independent, 
then the equation f(m, n) = k implies a finite number of systems of the form: 

aim + bin = di 
a 2 m + b 2 n = d 2 

where d\, d 2 run over the divisors^ of k. Each of these systems has at most one solution, 
implying that the equation f(m, n) = k has only a finite number of solutions. For the 
rest we assume that f(m,n) has no two linearly independent linear factors. Then the 
polynomial f(m, n) is not of the form c • g(m, n) 2 only in the following four cases: 

1. f(m, n) is irreducible 

2. f(m,m) = g 2 (m,n)h 2 (m,n) 

3. f(m,m) = gi{m,n)h^(m,n) 

4. f(m,m) = gi(m,n) 2 h 2 (m,n) 

where gi(m, n) and hi(m, n) are irreducible polynomials of degree i with integer coefficients; 
and g 2 (m, n) and h 2 (m, n) are linearly independent. 



4 Unless specified otherwise, the divisors of an integer include both positive and negative divisors. 



In the case 1, the polynomial f(m,n) is irreducible, implying that the Diophantine 
equation f(m, n) = k has a finite number of solutions due to the Thue theorem |14|. 

In the case 2, the equation f(m,n) = k is equivalent to the union of the following 
systems of equations: 

giim, n) = d 
h 2 (m,n) = k/d 

where d goes over the divisors of k. Each of these systems has only a finite number of 
solutions by Theorem |2] Therefore, the total number of solutions is also finite in this case. 

In the case 3, the solutions to f(tn, n) = k must also satisfy a Thue equation h 3 (m, n) = d 
for some divisor d of k. Since for each such d there is only a finite number of solutions, the 
total number of solutions to f(m, n) = k is finite. 

In the case 4, we let gikm, n) = g\(m, n) 2 and proceed as in the case 2 above. □ 

Theorem 4. A Diophantine equation f(m, n) = k, where k ^ and 

f(m,n) = A(D 2 m 4 + E 2 n A ) + 2Bmn(Dm 2 + En 2 ) + Cm 2 n 2 

with integer coefficients A, B, C, D, E, may have an infinite number of integer solutions only if 
1A 2 DE + B 2 - AC = 0. 

Proof According to Theorem [3J, the equation f(m, n) = k may have an infinite number of 
solutions only if 

f(m, n) = c ■ (urn 2 + vmn + wn 2 ) 2 = c(u 2 m 4 + w 2 n 4 ) + 2cvmn(um + wn) + c(2uw + v 2 )m 2 n 2 

for some integer coefficients c, u, v, w. In this case A = cr 2 , B = cvr, C = c(2uw+v 2 ), D = u/r, 
E = w/r for some non-zero rational number r and thus 

2A 2 DE + B 2 -AC = 2(cr 2 ) 2 (u/r)(w/r) + (cvr) 2 - (cr 2 )c(v 2 + 2uw) 
= c 2 r 2 (2uw + v 2 - (2uw + v 2 )) = 0. 

□ 

The theorem below represents a correct version of Corollary 6.3.8 in [2] which appears 
seriously flawedfl 

Theorem 5. Let A, B, C be non-zero integers and let (xo, yo, Zo) be a particular non-trivial integer 
solution of the Diophantine equation Ax 2 + By 2 + Cz 2 = with Zq + 0. The general integer 
solution to the above equation is given by 

V 

(x,y,z) = - (P x (m,n), P y (m,n), P z (m,n)) (3) 



5 Applying Corollary 6.3.8 [2| to the equation -25x 2 + y 1 = -z 2 with a particular solution (xo/i/0/ z o) 
( 1 , 0, 5), all integer solutions to this equation are given by (x, y, z) = \ (-25s 2 - t 2 , -50s£, -125s 2 + 5f 2 ), where 
the integers s, t as well as u, v are coprime and v \ 50. However, it is easy seen that, under the last restriction, 
the solution (x, y, z) = (1, 3, 4) cannot be obtained. 



where m, n as well as p, q are coprime integers with q > dividing 2 lcm(A, B)Cz 2 , and 

P x (m, n) = x Am 2 + 2y Bmn - x Bn 2 , 
P y {m, n) = -y Am 2 + 2x^Amn + y^Bn 2 , 
P z (m,n) = ZoAm 2 + zoBn 2 . 

Proof. The general solution © to the equation Ax 2 + By 2 = Cz 2 is given in Corollary 6.3.6 
in (2j. It remains to prove that q divides 2 lcm(A, B)Cz^. 

The values of x, y, and z in ^ are integer if and only if q divides each of P x (m,n), 
P y {m, n), and P z (m, n). Then q also divides 

A(y m - 2x n) ■ P x (m, n) + (x Am + 2y Bn) ■ P y (m, n) = 2BCz\k' 

as well as 

(2x Am + yoBn) ■ P x (m,n) + (-2y Bm + x Bn) ■ P y (m,n) = 2ACz\nv > . 
Therefore, q divides 

gcd(2BCz^ 3 ,2ACz> 3 ) = 2gcd(A / B)Czggcd( gcd( ^ g) n 3 / gcd ^ /B) ^ 3 ) 

which in turn divides 

□ 

Theorem 6. A system of Diophantine equations 

a\x 2 + b\y 2 + C\Z 2 = d\ 
a 2 x 2 + b 2 y 2 + c 2 z 2 = d 2 

where ai, hi, C{, di (i = 1, 2) are integers and in the matrix of coefficients 

fli b\ C\ d\ 
a 2 b 2 c 2 d 2 

every 2x2 minor is non-zero, has a finite number of solutions. 

Proof. It is clear that (di, d 2 ) ^ (0, 0). Without loss of generality we assume that d\ + 0. 

If x = then the resulting system has a finite number of solutions due to Theorem |2l 
For the rest of the proof assume that x + 0. 

Multiplying the first equation of (ID by d 2 and subtracting the second equation multi- 
plied by d\, we get a homogeneous Diophantine equation: 



M = 



Ax 2 + By 2 + Cz 2 = 



(5) 



where A = a\d 2 - a 2 d\, B = b\d 2 - b 2 d lr and C = C\d 2 - c 2 d\ represent minors of the matrix 
M and thus are non-zero^ 

Let (x , y , Zq) be a particular solution to equation © (for general methods of finding, 
particular solutions to quadratic homogeneous equations in three variables, see 12]] j3 
with Zo + 0, i.e., 

Ax\ + By 2 + Cz 2 =0, z + (6) 

then by Theorem[5]the general solution to equation © is given by ©, where gcd(m, n) = 1, 
gcd(p,^) = 1, and q > divides 2\cm{A,B)Cz 2 . 

Plugging the general solution © into the first equation of (JJ), we have 

f(m,n)=d 1 ft) (7) 



where 



f(m, n) d = {aix 2 + b\y\ + c 1 zl)(A 2 m 4 + B 2 n 4 ) 
+ 4x i/o( fl iS - b\A){Am 2 - Bn 2 )mn 
+ (^ylB 2 - 2AB{a x x\ + b x y\ - c x z\) + ^b x x\A 2 )m 2 n 2 . 

We remark that equation © represents a collection of Diophantine equations indexed 
by coprime integers p and q ranging over the divisors of respectively ^d\ld' x and 
2 lcm(A, B)Cz 2 , where d' x is the square-free core 

Theorem II] implies that equation © may have infinitely many solutions only if 

- 2(aiXq + biyl + c^lfAB + (2x i/o(«iB - hA)) 2 

- (aiXq + hyl + c^iAaxylB 2 - 2AB(aiX 2 + b x y\ - CiZ 2 Q ) + Ab^A 2 ) 
= - 4(Ab!xl + Bbiyl + Bc^iMiX^ + Ba x y\ + Acyfy 

= - 4z t (Bci - Cbi)(Ad - Cat) 
= - 4zQd 2 (biC 2 - b 2 C\)(a\c 2 - a 2 c\) 

is zero. However, each factor in this product is non-zero. Therefore, equation © has only 
a finite number of integer solutions. □ 



6 From computational perspective, it is worth to divide each of the coefficients A, B, C by their greatest 
common divisor gcd(A, B, C) before solving equation 10. 

7 PARI/GP provides a function bnfisnorm, which can be used to compute a particular solution to a quadratic 
homogeneous equation in three indeterminants. 

8 Since the number of resulting equations is proportional to the number of divisors of 2 lcm(A, B)CZq, 
from computational perspective it may be worth to choose a particular solution (xo, yo, Zq) and re-label the 
variables (if needed) such that the number of divisors of 2 lcm(A, B)Cz;? is minimized. 



3 Intersection of two Lucas sequences 



3.1 Intersection of two IZ-sequences 

Ifanintegerxbelongstotheintersection{Lr )i (Pi,Qi)}and{LT )3 (P2/Q2)}/Where|Qi| = IQ2I = 1 
and (Pi, Qi) ^ (P2, Qi), then dTJ) implies that x satisfies a system of Diophantine equations: 

y 2 - D x x 2 = ±4 

[z 2 - D 2 x 2 = ±4 K ' 

where Di = P 2 - 4Qi and D 2 = P 2 - 4Q 2 . The matrix of its coefficients 

-D x 1 ±4 
-D 2 1+4 

has all 2 x 2 minors non-zero unless Di = D 2 and the signs of 4's are the same. 

If Di ^ D 2 , then system ((S) has only a finite number of solutions by Theorem [6j We 
note that it always has solutions (x,y,z) = (0,2,2) and (x,y,z) = (l,Pi,P 2 ) corresponding 
to common elements = U (P lf Qj) = U (P 2 , Q 2 ) and 1 = U x (P lf Q x ) = U X {P 2 , Q 2 ). 

The case of D\ — D 2 is addressed by the following theorem. 

Theorem 7. For distinct non-degenerate sequences {U n (Pi,Qi)} and {U„(P 2 , Q 2 )} wz'f/i P, > 
and \Qi\ = 1 (i = 1,2), the following conditions are equivalent: 

• their intersection is infinite; 

• their discriminants are equal; 



these sequences are Fibonacci numbers {U n (l, -1)} and their bisection {U n (3, 1)} 
{U 2n (l,-1)\. 



Proof As discussed above, the intersection of {LT„(Pi, Qi)} and {U n (P 2 , Q 2 )} can be infinite 
only if their discriminants are equal, i.e., P 2 - 4Qi = P\ - 4Q 2 . 

If Qi = Q 2 then Pi = P 2 , which is not the case for distinct sequences. 

If q 1 = -Q 2 then P 2 - P\ = ±8, implying that (P lf P 2 ) = (3, 1) or (P lf P 2 ) = (1, 3) and thus 
the sequences are {U n (l,-1)\ and {U n (3, 1)}. Since U n (3, 1) = U 2n (l,—1) for every integer 
n, {U n (3, 1)} represents a subsequence of U n (l,— 1), implying that they have an infinite 
intersection: 

{LT n (l / -l)}n{LI n (3 / l)} = {LJ„(3,l)}. 

□ 



3.2 Intersection of a IT-sequence and a V-sequence 

If an integer x belongs to the intersection {U n (Pi, Qi)} and {V n (P 2 , Q 2 )), where |Qi I = IQ2 
1, then ([I]) implies that x satisfies a system of Diophantine equations: 

y 2 - D x x 2 = ±4 
x 2 - D 2 z 2 = ±4 



where D\ = P? — 4Qi and D 2 — P 2 - 4Q 2 . The matrix of its coefficients 

"-Di 1 ±4 
1 -D 2 ±4 

has all 2 X 2 minors non-zero unless Di = 1, which is not possible. By Theorem [6l system 
^ has only a finite number of solutions. 



3.3 Intersection of two V-sequences 

Ifanintegerxbelongstotheintersection{ 1 l/ fI (P 1 ,Q 1 )}and{V r f ,(P2,Q2)},where|Qi| = \Q 2 \ = 1 
and (Pi, Qi) ^ (P 2 , Q 2 ), then §\§ implies that x satisfies a system of Diophantine equations: 

[j 2 - D 2 z 2 = ±4 1 ' 

where D\ = P 2 - 4Qi and D 2 = P 2 - 4Q 2 . The matrix of its coefficients 

1 -D a ±4 
1 -D 2 ±4 

has all 2 X 2 minors non-zero unless the signs of 4's in the r.h.s. are the same. In the 
case of distinct signs in the r.h.s. of system ((10) , it has a finite number of solutions by 
Theorem [6j Below we address the case of the same signs, which leads to the following 
system of equation: 

D,y 2 -D 2 z 2 = 0. ( ' 

The second equation of system ([iT) always has a solution y = z = corresponding 
to the common element 2 of any two ^/-sequences. On the other hand, the system has 
nonzero solutions only if the product DiD 2 is a square. 

Theorem 8. The intersection of distinct non-degenerate sequences V n (Pi, Qi) and V n {P 2 , Q 2 ) 
with Pi > and |Q,| = 1 (i = 1,2) is infinite if and only if the product of their discriminants is 
a square. Moreover, in this case the intersection itself represents a Lucas sequence V n (P',Q')for 
some P' > and \Q'\ = 1. 

Proof. As shown above, the intersection of V n (Pi, Qi) and V n (P 2 , Q 2 ) can be infinite only if 
the product of their discriminants D\ = P 7 : - 4Q X and D 2 = P\ - 4Q 2 is a square. To prove 
the converse, assume that D\D 2 is a square, i.e., D\ = c&\ and D 2 = cd\ for some integers 
d\, d 2 with gcd(rfi, d 2 ) = 1, and non-square integer c. 



We first note that in this case, system (|10|) with opposite signs in the r.h.s. does not 
have a solution. Indeed, it implies that 



x 4 - 16 = (x 2 + 4)(x 2 - 4) = D^yh 2 = {cd x d 2 yzf 



and hence 

(x 2 - cd\d 2 yz){x 2 + cd\d 2 yz) = 16 



whose the only solutions x = ±2 do not satisfy the original system ((IP)) . Therefore, the 
intersection of {V n (Pi, Qi)} and {V n (P 2 , Q 2 )} necessarily represents the solutions to ((10) with 
the same signs in the r.h.s., which is equivalent to (|TT|) . 



The general solution to the second equation of (|Tlj) is (y,z) = (±d 2 t,±d 1 t) for £ e Z, 
where the signs are independent of each other. Plugging it into the first equation of (|TT]) , 
we get a Pellian equation: 



x 2 - Df 2 = ±4 (12) 

Where D = gcdXoa) = Cd ? rf 2" 

Notice that both D\ and D 2 are congruent to or 1 modulo 4, and so is D. Therefore, it 
has infinitely many solutions (x, t) = (Xk, t k ) of the form: 

x k + t k Vd (xi + ti Vd\' 1 



ke Z, 
z \ z / 

where (xi, fi) is a solution with the smallest positive h (see Q, Proposition 6.3.16). It 
follows that the sequence x k is simply Vk(h, sign(x 2 - Df 2 )). 

Therefore, the intersection of Lucas sequences V n (Pi,Qi) and V n (P 2 , Q 2 ) forms Lucas 
sequence V k (ti, sign(x 2 - Dt 2 )). 

□ 



4 Computing intersections 

We note that the theorems in Section [T] have constructive proofs that allow straightfor- 
ward implementation. A bottleneck of our approach in practice is Thue equations that 
are somewhat hard to solve. Nevertheless, there exist algorithms for solving Thue equa- 
tions IHBIHI, while PARI/GP computer algebra system ||T3l provides an implementation of 
such an algorithm. Below we present some numerical results obtained with the help of 
PARI/GP. 

4.1 Intersection of Fibonacci and Pell numbers 

Theorem 9. The only numbers that are both Fibonacci and Pell are 0, 1, 2, and 5. 

Proof. To find the intersection of Fibonacci numbers li n (l, — 1) and Pell numbers IT„(2, -1), 
we first compute their discriminants Di = 5 and D 2 = 8 and plug them into system ©: 

V-5x 2 = ±4 

|z 2 -8x 2 = ±4 {L) 



As it was shown in Section 13.11 this system has only a finite number of solutions. We 
will find them all, following the outline of the proof of Theorem [6j First we obtain a 



homogeneous equation (depending on the choice of signs of 4's in the r.h.s.), resulting in 
two possible systems of equations: 



y 2 - 5x 2 = ±4 
\3x 2 + y 2 -z 2 = 



y 2 - 5x 2 = ±4 

-13x 2 + y 2 + z 2 = 



that we solve below. 

Using a particular solution (xq, yo,Zo) = (0, 1, 1) to the equation 3x 2 + y 2 - z 2 = 0, from 
© we obtain a general solution: 

! x = 2mn- 
3/ = {3m 2 -n 2 ) p - 
z = (3m 2 + n 2 ) p - 



where the integers m, n as well as p, q are coprime, and q > divides 6. Plugging these 

jua 



expressions into the first equation of (fl3|) , we obtain equations 

9/;/ 4 -26//r/r + /r = ±\ — 

V 



where p divides 2. Since the polynomial in the l.h.s. is irreducible, these equations are 
Thue equations. 

According to PARI/GP, they have the following coprime solutions (m,n) and corre- 
sponding triples (x, y, z) (up to signs): 





(M,M) 


(W/lyLN) 


(2,1) 


(0,1) 


(0,2,2) 


(1,2) 


(0,2), (1,1), (1,5) 


(0,2,2), (1,1,2), (5,11,14) 


(2,3) 


(1,0) 


(0,2,2) 


(1,6) 


(2,0), (1,3) 


(0,2,2), (1,1,2), (5,11,14) 



Using a particular solution (xo, yo, Zq) — (1,2, 3) to the equation -13x 2 + y 2 + z 2 = 0, from 
© we obtain a general solution: 

x = (13m 2 - 4mn + n 2 )^ 
y = (26m 2 - 26mn + 2n 2 ) v - 
[z = (39m 2 - 3n 2 ) v - 

where the integers m, n as well as p, q are coprime, and q > divides 234. Plugging these 
expressions into the first equation of (fi~3~l) , we obtain equations 

-169m 4 - 832ra 3 n + 570m 2 n 2 - 64mn 3 - n 4 = + |y J 

where p divides 2. Since the polynomial in the l.h.s. is irreducible, these equations are 
Thue equations. According to PARI/GP, they have the following coprime solutions (m,n) 
(up to simultaneous change of signs) and corresponding triples (\x\, \y\, \z\): 



/ \ 

(p,q) 


(m,n) 


/III I I l\ 

(\x\,\y\, \z\) 


(2,1) 


(0,1) 


(2,4,6) 


(1/2) 


(0,2) 


(2,4,6) 


(2,9) 


(0,3), (1,2) 


(2,4,6) 


(2,13) 


(1,0) 


(2,4,6) 


(1,18) 


(0,6), (1,-1), (1,5), (2,4) 


(1,3,2), (2,4,6) 


(1,26) 


(2,0) 


(2,4,6) 


(2,117) 


(2,13), (3,0) 


(2,4,6) 


(1,234) 


(1,-13), (5, 13), (4, 26), (6,0) 


(1,3,2), (2,4,6) 



□ 

4.2 Intersections of some Lucas sequences 

In this section we lists some interesting numerical results and pose related open questions. 
As before, we consider only on non-degenerate sequence and in addition exclude the 
sequences with (P, Q) = (3, 1) corresponding to bisections of those with (P, Q) = (1, -1). 

We call the intersection of sequences 14 (Pi, Q\) and 14 (P2, Q2) non-trivial if it contains 
an integer m > 1 different from Pi and P2. Examples of non- trivial finite intersections are 

{14(3,-1)} n {14(19,1)} = {0,1,360}, 
{14(8,-1)} n {14(23, 1)} = {0,1,528}, 
{14(3,-1)} n {14(6,1)} = {0,1,1189}. 

For an integer m > 1, we can find all values of P > such that m e {U n (P, Q)\ for some 
Q = ±1 by solving the equation (followed from Theorem [T]) 

y 2 - m 2 p 2 = ±4:(m 2 ± 1). 

In particular, the only integers m < 10 9 that belong to two distinct sequences U. n {P\, Qi) 
and 14(P 2 , Q2) for some Pi ±m,P 1 + m, and |Qi| = |Q 2 | = 1, are 5, 360, 528, and 1189. Are 
there any integers m > 10 9 with the analogous property? 
Another notable intersection is 

{14(1, -l)}n {14(2,-1)} = {0,1,2,5} 

which is the only intersection of size four that we are aware about. Does there exist any 
other? 

Similarly, we call the intersection of U n (Pi, Qi) and V n (P2, Q2) non-trivial if it contains 
an integer m > 2 different from Pi, and P 2 . Examples of non-trivial intersections are 

{14(1, -l)} n {14(2, -l)} = {14(1, -l)} n {V n (6, l)} = {2, 34}, 

{U n (9,-l)}n{y„(2,-l)} = {82}. 

In particular, the only integers m, 2 < m < 10 9 , that belong the intersection U n {P\, Q\) and 
V n (P2, Qi) for some P 1 +m,P 2 + m, and |Qi| = |Q 2 | = 1, are 3, 29, 34, and 82. Are there any 
integers m > 10 9 with the analogous property? 



The intersections 

{lT B (l,-l)}n{V B (l,-l)} = {1,2,3}, 

{LT n (2,-l)}n{y„(l,-l)} = {1,2,29} 

have size 3. Are there any other? 

Again, we call the intersection of V n (P\, Qi) and V„(P 2 , Qi) non-trivial if it contains an 
integer different from m > 2, Pi, and P2. Examples of non-trivial finite intersections are 

{V„(l,-l)}n{V K (3,-l)} = {2,3,11}, 

{v„(2,-i)}n{v n (4,i)} = {2,i4}. 

Are there any other non-trivial finite intersections? 

Tablefflgives all parameters (Pi,Qi) and (P 2/ Q 2 ) such that Pi < P 2 < 1000, |Qi| = |Q 2 | = 1 
and {V n (P lf Qi)} n {V n (P 2 , Q 2 )\ is infinite. 

5 Lucas sequences with arbitrary initial terms 

For fixed P, Q e Z, |Q| = 1, consider a sequence that satisfy the same recurrent relation as 
U„(P, Q) but with initial terms a, b £ Z: 

Uf(P,Q) = a, Uf(P,Q) = b, U a *JP,Q) = P-U%\P,Q)-Q-Uf_ 1 (P,Q),n>l. 

It is easy to see that it can be expressed in terms of conventional Lucas sequences as 
follows: 

U a /(P, Q) = (b- aP) ■ U n (P, Q) + a- U n+1 (P, Q). 
This representation allows us to define a counterpart sequence: 

V a M Q) "= (b - aP) ■ V n {P, Q) + a- V n+l {P, Q) = U^^iP, Q). 

It is easy to see that the generalized Lucas sequences satisfy an identity similar to dU): 

Vf (P, Q) 2 - D • Uf{P t Qf = ((2b - aP) 2 - Da 2 )Q n 

where D = P 2 - 4Q. However, there is no general analog of Theorem [T] and hence the 
identity 

y 2 -D-x 2 = ±((2b - aP) 1 - Da 2 ) (14) 

can be used only as a necessary condition for testing membership of x in W t ; h (P l Q). In 
particular, while (x, y) = (8, 34) is a solution to equation (p~4|> for (P, Q, a, b) = (1, -1, 1, 15), 
8 is not an element of U); 15 (l, -1). 

The quadratic equation <[T4|) allows one to compute the intersection of Lucas sequences 
with arbitrary initial terms and prove its finiteness similarly to the case of classic Lucas 
sequences. 
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(123,1) 


(4,-1) 
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(76,-1) 
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(141422324,-1) 




(12,1) 
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(142,1) 


(4,-1) 


(843,1) 


(599074578,1) 




(12,-1) 


(146,1) 


(146,1) 



Table 1: Parameters such that {V n (P lf Qi)} n {V n (P 2 , Q 2 )} = {V„(P, Q)}. 
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(167 11 
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Table 1: (con'd) Parameters such that {V n (P lf Q 1 )} n {V n (P 2 , Q 2 )} = {V n (P, Q)}. 



